www.gusucode.com > 基于matlab编程蒙特卡罗仿真设计源码程序 > 基于matlab编程蒙特卡罗仿真设计源码程序/第五讲 蒙特卡洛方法的应用/sjtdf_pi2.m
function piguji=sjtdf_pi2(mm) %mm 是随机实验次数 frq=0; xrandnum = unifrnd(0,1,1,mm); yrandnum = unifrnd(0,1,1,mm); for ii=1:mm if (sqrt(1-xrandnum(1,ii)^2)>=yrandnum(1,ii)) frq=frq+1; end end piguji=4*frq/mm